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Abstract 

We investigate the evolution of a single qubit subject to a contin- 
uous unitary dynamics and an additional interrupting influence which 
occurs periodically. One may imagine a dynamically evolving closed 
quantum system which becomes open at certain times. The inter- 
rupting influence is represented by an operation, which is assumed to 
equivalently describe a non-selective unsharp measurement. It may be 
decomposed into a positive operator, which in case of a measurement 
represents the pure measurement part, followed by an unitary back- 
action operator. Equations of motion for the state evolution are derived 
in the form of difference equations. It is shown that the "free" Hamil- 
tonian is completed by an averaged Hamiltonian, which goes back to 
the back-action. The positive operator specifies a decoherence rate and 
results in a decoherence term. The continuum limit to a master equa- 
tion is performed. The selective evolution is discussed and correcting 
higher order terms are worked out in an Appendix. 



1 Introduction 

Experimental and theoretical studies of the dynamics of single two-level 
systems have become very important in the context of quantum computation 
and quantum information. In this article we investigate the evolution of a 
single qubit subject to a continuously acting unitary dynamics (undisturbed 
or "free" dynamics) with operator U as well as affected by an interrupting 
additional influence, which is non-unitary and acts periodically at times 
tn = to + nr, n = 1 , 2 .... The duration 6t of this influence is assumed to 
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be much shorter than r, so that it can be neglected. One may imagine a 
dynamically evolving closed quantum system which becomes open at times 
tn- The corresponding single influence is represented by an operation £ 
which transforms the state of the qubit given by its reduced density operator 
p according to 

p ^ Sip) = MkpM+ (1) 
k=± 

with operation elements Mj, , which are sometimes also called Kraus oper- 
ators. £ is assumed to be trace preserving. The representation (||) of the 
single operation is called the operator-sum representation or Kraus repre- 
sentation. 

Such a periodically occurring, nearly instantaneous change can be caused 
by a recurring interaction with a second system provided that this system 
does not "remember" the influence it may have experienced from the qubit 
at former times (Markov process). Typically the second system could be 
an environment or it could consist of a number of systems of the same kind 
which interact only once with the qubit, as it is the case in a sequence of 
scattering processes. 

According to eqn. (|l|) we are not dealing with the most general form of 
such an influence on a qubit, which would correspond to an operation with 
four operation elements, but we restrict to interactions which may be repre- 
sented by only two elements M_|_ and M_. Pairs of operation elements can 
describe such important operations as for example amplitude damping and 
phase damping, bit flips and phase flips as well as projection measurements 
and unsharp measurements. The concept of an unsharp measurement will be 
explained below. We will restrict to an operation which is equivalent, as far 
as the map £ is concerned, to a non-selective unsharp measurement. Note 
that this may still comprise many physical processes which at flrst glance 
do not look like an unsharp measurement. An example is the experiment 
of Brune, Haroche et al. to measure the number of photons in a cavity, 
which is discussed under the aspect of an unsharp measurement in |2| . 

According to the polar decomposition theorem each operation element 
M± may be written as a product of a unitary operator and a positive oper- 
ator 

M± = C/±|M±|. (2) 

We are dealing with a class of generalized measurements with two outcomes 
-|- and — . The unitary part U± extracts no information from the qubit. It is 
in a different context often called the feedback part of the quantum operation 
[P, ^. The set {|M_p, |M_|_p} represents a positive-operator valued measure 
(POVM). The demand of unsharpness of the measurement is introduced 
below in specifying these \M±\ further. A sequence of such generalized 
measurements are of practical importance because they can be employed to 
explore the original dynamics of the system |5| or to control its dynamics 
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by means of a specific feedback 

In the non-selective case of eqn. (|l|), wlien tlie measurement results it are 
not read off, our total physical setup defined by U, r, U± and M± may also 
be regarded as a particular noisy channel. Below the unitary operators U± 
will not be neglected because they form an important part of the sequential 
operations in realistic situations. In general it would need a non-trivial 
feedback procedure to eliminate their influence, compare for an example. 

It is our goal to derive equations of motions for the state of the system 
subject to the periodical influence. We will do this in form of difference 
equations and - on a coarse grained time scale - in form of master equations. 
Difference equations take into account the discrete nature of the influence 
due to the finite time r between its occurrences and grant therefore a more 
exact description. It may reveal specific physical influences which cannot be 
seen in the continuum limit furthermore. 

We also look at the limit r — > of continuous measurements and com- 
pute the corresponding master equation (non-selective description). The 
master equation can be understood as a means to compute approximately 
the dynamics of a system which is subject to a sequence of operations of the 
type (|^). Since there are a lot more mathematical methods to solve differ- 
ential equations than difference equations, it is often useful to approximate 
sequential measurements by continuous measurements. 

Master equations for special cases of measurements with non-minimal 
disturbance of the state , i.e. with a non-vanishing unitary part of the op- 
eration elements have been considered in the literature. A master equation 
for general feedback was derived by Wiseman Q. In the Markovian limit, 
if the time delay between feedback and measurement vanishes and the feed- 
back depends only on the outcome of the last measurement (instantaneous 
feedback), the action of the feedback can be represented as unitary part 
of the operation of the measurement. However Q does not comprise our 
results since it deals with a special kind of continuous measurements. They 
have poissonian statistics and allow finite state changes during infinitesimal 
time intervals. More precisely this means that only very seldom a certain 
measurement result occurs which is then connected to a finite state change 
during a infinitesimal time while for other measurement results the state 
changes only infinitesimally. We are excluding Poissonian statistics and re- 
quire the state change in the continuum limit to be infinitesimal during 
infinitesimal times. Thus Wiseman's and our studies do not overlap. 

In a later paper Wiseman |Q] employed the operation formalism to an- 
alyze a homodyne measurement in quantum optics to apply instantaneous 
feedback in order to minimize disturbance, i.e., compensate the unitary part 
of the operation. Korotkov investigated a measurement with non-minimal 
disturbance in the context of continuous measurement of a qubit by means 
of a single electron transistor ||8|, ^ . He noted that this non-minimal distur- 
bance acts in the master equation like a change of the distance between the 
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two energy levels of the qubit. We find this effect as a special case of our 
studies {[U±,H] = 0). Korotkov also derives a modification of the stochastic 
master equation (selective regime) due to the non- minimal disturbance 
which is averaged out in the non-selective case. 

Many examples of the application of instantaneous feedback in contin- 
uous measurements by means of external changes of the Hamiltonian of 
the system can be found in the literature of quantum control (e.g. ^). 
They correspond to special choices of the unitary part of the measurement 
operation and obey the equations of motions derived here, provided the 
measurements do not inflict finite changes during infinitesimal times. 

In the context of quantum dissipation the influence of a heath bath on a 
infinite dimensional quantum system has been investigated by Caldeira and 
Leggett . The coupling was such that in terms of the operation formalism 
the corresponding operation of the system had a unitary part additional to 
the one stemming form its free original evolution. They derived a master 
equation for high temperatures which was later modified to also describe 
medium temperatures |11|. Although we are looking at a qubit we find in 
the difference equation among others similar terms like there, but only one 
of them survives in the continuum limit. 

In order to discuss the problem in an illustrative way but without re- 
striction of generality, we treat the periodical influence in terms of a se- 
quential measurement. We proceed as follows. We first consider the sin- 
gle measurements of the sequence. Then we bundle the whole sequence in 
subsequences of N measurements ("N-series"). The resulting operation has 
Gaussian shape. Afterwards we integrate over all outcomes to obtain the av- 
erage state change due to a N-series (non-selective regime) . We then discuss 
how to find the right continuum limit which conserves the physical charac- 
teristics of the N-series and derive the master equation. Finally we deal with 
the selective regime of measurement and write down the stochastic master 
equation. The Appendices serve to derive the difference equation for the 
non-selective regime up to second order in the occurring small parameters. 



2 The single quantum operation 

An example for a physical realisation of the operation £ from (||) is given 
by a qubit which interacts at times t„ unitarily for a short duration 6t with 
an environment and thus becomes an open system. The resulting change 
of its reduced density operator may formally be expressed with the help of 
operation elements M± as in eqn. (HI). Just the same change of the reduced 
density operator results if i) a projection measurement on the environment 
is performed with outcome -|- or — transferring the qubit to the states 

p ^ p^ = M±pM± . (3) 
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respectively and ii) the outcome is not read off (non-selective case). In 
the generic case, eqn. (^) describes thereby a generalized measurement of 
the qubit. For simplicity reasons the terminology we are going to use will 
refer to measurements, but the results apply equally to any operation with 
operator-sum representation (0) if the same specifications of M± are made. 
This is independent of how the operation is experimentally realized. 

Because of the polar decomposition theorem, the operation elements M± 
may be written as products of a unitary operator and a positive operator 

M± = C/±|M±|. (4) 

We introduce the POVM effects 

E± = \M±\\ (5) 

which obey the completeness relation 

E+ + E^ = l. (6) 

The probability of the outcome + or — is given by 

P± = {E±)p (7) 

with |)_|_ + p_ = 1. 

Eqn. represents a decomposition of the operation into a pure mea- 
surement part described by |M-|-|, followed by a unitary back-action given 
by U± depending on the result + or — . These denominations are justi- 
fied for the following reasons: All the information which can be read off 
from the meter is related to \M±\ which therefore represents the unavoid- 
able minimal disturbance. The unitary operators leave the von Neumann 
entropy unchanged and therefore do not allow to export information to an 
observer. Because they depend on the result + or — , they may be inter- 
preted as a specific back-action of the measuring apparatus inducing an 
additional Hamiltonian evolution of the qubit. We formally introduce the 
corresponding Hamiltonians -ff-j- according to 

U± =: exp {-^H±T^ . (8) 

This unitary back-action represents an important part of the quantum op- 
eration and appears naturally in the generic situation. 

According to eqns. (|5|) and (^) |M+p and |Af_p commute. Therefore 
we can find orthonormal basisvectors |1) and |2) of the qubit Hilbertspace 
with respect to which |M±p are diagonal. We introduce the eigenvalues pi 
and p2 of |M+p, which are positive and because of eqn. obey < pi^2 < 
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1. Without restriction of generahty we choose p2 > pi- Reading off the 
eigenvalues of |M_p from (^) and taking the square root we find 

\M+\ := v ^|l)(l | + v^|2 )(2| (9) 
|M_| := yr^|l)(l| + yT^|2)(2| 

The elements |M_|_| and |M_| commute. We will characterize the operation 
later on by the parameters 

Po ■■= ^{Pi + P2) , Ap := p2 - pi (10) 
with < Ap < 1. Introducing 

a,:=|l)(l|-|2)(2| (11) 
the effects E± of eqn. (|5|) are rewritten in the form 

E+ = Pol - ^Apa„ E_ = il- po)l + ^Apa, . (12) 

In the limiting case Ap = 1, the pure part of the measurement (^) results 
in a projection on |1) or |2) depending on the measurement outcome + or — . 
We call this a sharp measurement of an observable with eigenvectors |1) and 
1 2), for example az- Note that also for a sharp measurement the result of 
the quantum operation (jl]) will in general not be the state |1) or |2) because 
of the remaining influence of the unitary back-action. 

In the contrary limit Ap <C 1 the |M_|_| and |-/Vf_| are nearly propor- 
tional to the identity operator. The probability p-^- (or p_) to obtain the 
measurement result + (or — ) is then nearly independent of the initial state 
of the qubit. There is almost no state discrimination. Because of this low 
sensitivity we call this an unsharp measurement. Note that in this limit the 
parameters po and I — po become, because of eqns. and (12), approx- 
imately the mean probabilities to obtain the measurement results + or — 
respectively. "Unsharpness" does of course not originate from a measure- 
ment apparatus which is "broken". 

We are now able to further specify the particular class of quantum op- 
erations (|l]) which we are going to discuss. We will restrict to the case that 
the pure measurement part represents an unsharp measurement: Ap <C 1. 

We will mainly be interested in the non-selectiv case where an informa- 
tion about the results it or the corresponding states of the environment is 
not available. The influence on the qubit at times t„ may then be written 
in the operator-sum representation as in eqn. (|l]). 

p^£{p)=J2 MkpM+, (13) 
k=± 

whereby 

J2M+Mk = l. (14) 

k=± 

because of (g). The quantum operations at times tn are trace preserving. 
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3 N-series and related operation 



The time between measurements is r. We assume that the duration 5t of a 
measurement is much shorter than r. The undisturbed or "free" dynamics 
of the system between the measurements is given by the Hamiltonian H. We 
bundle N consecutive measurements to a N-series of duration At = Nt as 



we have done in [Q] (cp. also [|12[)- This procedure has several advantages. 
We will obtain effects of a Gaussian structure. This enables us to work out 
the operator sum explicity. A comparison with the results in the literature 
regarding continuous measurements becomes more evident. And finally the 
discussion of the selective case is simpler. 
We require 

iV>l. (15) 
We relate N to the sharpness Ap ^ 1 of the measurement by 

N ■Ap<t.l (16) 



and demand in addition 
and 



At\\H\\<^h (17) 



At||iJ±|| < (18) 

This means that the influence of the undisturbed dynamics of the qubit and 
the unitary back-action dynamics due to the measurements are both small 
over the duration At of a N-series. With At = Nt we have obtained above 
restrictions for A^, r, H and H±. 

The density operator resulting at the end of a N-series of measurements 
with results mi, ...rriN, each of which can assume the values "+" and "— ", 
read 

pit + At) = Mm^U...Mrn,Up{t)U+M+^ ... U+ (19) 

with 

U :=exp{-^HT} . (20) 

The influences of the systems dynamics and the measurement will in 
general not commute. The following relation is derived in Appendix A: 

M^^UMm,,_,U ...Mm,U = M^^Mm^_, ...Mrn,U^{l + Cl) (21) 

with 

IICill < 0{NApAt\\H\\/n) + 0{At^\\H\\ ui&x{\\H±\\} /Ti^) . (22) 

Also the operations M± will not commute. Based on the decomposition (^) 
we show in Appendix A that 

C/^^IM^^I ... [/^JM^JC/^ = C/i^+[/^-^+|M+r+|M_r-^+^^(l + C2) . 

(23) 
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with 



IC2II < 0{N /\pM\\m.&-^{\\H±\\} /h) + 0{I\t^\\H+\\\\H_\\/n^ 



(24) 



and — are the total numbers of measurement results "+" and "— " 
in the N-series respectively. Because of the assumptions (|l^), ( [l7| ) and (18) 
we may neglect Ci and C2. The calculation which takes into account C\ 
and C2, is given in the Appendix C. 

In our approximation the influence ( |^ of the N-series becomes a func- 
tion of Nj^ only, independent of the (^^) different orderings of the "+" 
and "— " results (cp. eqn. (p^)). Therefore the total N-series of duration 
Ai including the "free" dynamics corresponds to a quantum operation with 
operation elements 



M{N^,N) = U^+U^-^+\M{N+,N)\ exp { - ^H^t] 



(25) 



with 

\M{N+,N)\ 



' N 



|M^r+|M_r-^+ 



(26) 



N+ JV-JV, 
[Pl' (1-Pl)^|l)(l|+P2' (I-P2) 



|2)(2| 



The square root in front ensures the completeness relation of the effects 

AT 

M{N+, N)+M{N+, N) = l. (27) 

N+=0 

The unitary back-action part caused by the N-series measurements can be 
expressed as 



U{N+,N) 



exp |-^(7V+F+ + {N- N+)H^)t^ . 



(28) 



We now make use of the condition that N is large, so that \M(N^, N)\ 
of eqn. (pi) may approximately be written in form of a Gaussian 



\M{N+,N)\ 



^2tiNE+E^ 



exp 



(29) 



which contains the operators of (12). Because we assumed that the mea- 
surements are unsharp and therefore Ap <^ 1, the "spread" of the Gaussian 
becomes in lowest order a c-number 

E+E_ = {poil - po) - \i^pf) 1 + Ap(po - b'^z = Po(l - Po) , (30) 
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where we have ignored terms of order Ap and higher on the right-hand side 
of (30). The error thus committed in the Gaussian in (^) is of order NAp^, 
which can be seen by inserting from ( [T^ ) and N^/N from (^Tj) and 
expanding Gaussian (^) in powers of Ap. A more detailed calculation can 
be found in Appendix C. 

We introduce a new variable s to replace the readout N^/N according 

to 

N, 1 

^ =: po - ^^ps . (31) 

Because is large we may approximately regard s to be continuous. Its 
range is limited by: 

0<po- ^Aps < 1 . (32) 
In addition we introduce the new quantity 

(Apf 

7 := -j — ir. ^ , (33) 

which will turn out to be the decoherence rate. It contains apart from Ap 
and pq also the time interval r between two measurements. These three 
parameters characterize completely the influence of the sequence of pure 
measurements. 7 increases when the measurements become sharper and 
accordingly have a stronger influence on the qubit. A decreasing time differ- 
ence T between two measurements results as well in an increase of 7. This 
reflects a Zeno type effect which also happens for unsharp measurements. 

Installing 7 we get the ultimate form of the operations, valid for Ap <C 
T/At < 1: 

1 ( {C^z- S)^ 



IMJ = , exp <^ - 7 ^ \ ' At \ . (34) 

' ^27r/(7At) ^1^4 j ^ ^ 

The resulting effects Eg = |Msp have Gaussian form. They show the char- 
acteristics which are known for instance from an unsharp position measure- 
ment as investigated e.g. in Instead of a continuous observable however, 
we are dealing here with a discrete observable. 

With reference to s, the complete operation element including the back- 
action and the "free" evolution is given by 

Ms = Us\Ms\ exp{-^i7At} , (35) 

where, using (|28|) and (|3l|), we obtain for Us which replaces U{Nj^,N): 

Us = exp j-^^AvAt - ^AHsApAt^ . (36) 

We have thereby introduced the averaged back-action Hamiltonian Hav and 
the difference AH of the back-action Hamiltonians, respectively: 

Hav ■■= PoH+ + {1 - po)H_ (37) 
AH := {H_ - H+) . (38) 
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Ms of ( |35| ) replaces M{N^,N) of eqn. (^) for the continuous variable s. 

The N-series operation elements above correspond to a continuous set of 
effects with the Gaussian distribution function 



Ps = {M+Ms)p . (39) 

The completeness relation is satisfied if we extend the range of s to the whole 
real axis: 

/oo 
M]M, ds = l. (40) 
-oo 

The statistical weight of the unphysical values of s will be negligible provided 
that po 7^ , 1 and this justifies the formal extension of the values of s beyond 
their physical range (p^). 



4 Non-selective evolution 

In the non-selective case the state change during a N-series can be expressed 
in the operator-sum representation as 

/oo 
MspM+ds . (41) 
-oo 

We are going to expand the r.h.s. up to linear terms in At. 

The unitary parts of the operation Ms which are generated by H , Hav 
and AH lead to 

Ap := p{t + At) - pit) = -'-[H + Hav, p{t) ]At + D{p{t)) - p{t) . (42) 



ds 



The integral 

D{p):= J exp^^-^AHsApAt^\Ms\p\Ms\ expi^^AHsApAt^ 

(43) 

over the parts which depend on s will first be calculated and then expanded. 

Introducing operators which act from the left and are denoted with L 
(e.g. a^p := (7^ p) as well as operators which act from the right and are 
denoted with R (e.g. cr^p := pcTz), we rewrite the integrand of D{p) 

xexp{-J((<Tj'-s)2 + (a?--»)2) Atj dsp. (44) 

It is important to take into account a further operator ordering for the 
integrand, namely that AH^,AH^ should remain leftmost and rightmost, 
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respectively. The resulting integral is Gaussian in s. It may be solved in a 
closed form: 



D{p) = exp 



A( p. (45) 



Now we expand it up to the leading linear term in At and restore the 
usual operator formalism according to the rules for L and R. This leads for 
example to 

= ia^-a^)[a,,p] = [a,,[a,,p]] (46) 

and all together to 

Dip) -p= (47) 



lWzA^z,p]]- ^ [AH, [AH, p]]-i^ [AH, {a,,p}]]At. 



While the second and the third term are proportional to small quantities 
(cp. (|l^) - ( [TsD) , the first contains the ratio of the two small quantities Ap^ 
and r (cp. (p3|)). We assume 7 not to be small. We thus obtain as final 
result for the state evolution during one N-series 

Ap=[-^[H + HAv,p]-l[^.,[<^z,p]])At. (48) 

The first term on the r.h.s. represents the unitary dynamical evolution re- 
lated to the "free" Hamiltonian H and to the averaged Hamiltonian H^y = 
PqH^ + (1 — Po)H^. po and 1 — po are approximately the probabilities that 
the back-action causes a Hamiltonian development with or H^ respec- 
tively (cp. section |2|). The second term on the r.h.s. reflects the decoherence 
induced by the pure measurement part [M±\ of the operation. The structure 
of both terms is clearly what one would expect on physical grounds. 

The second and third term of eqn. (^^ indicate additional physical ef- 
fects, which are to be expected in a higher order approximation. The second 
term corresponds to further decoherence induced by unitary back-action. 
The third term goes back to friction caused also by the back-action. A 
complete list of additional terms is given in Appendix C 

Eqn. (^) has been derived on the basis of the following approximations: 
We have changed the order of U, U± and [M±\ in the operation of the N- 
series, neglecting the commutators between them, cp. ( pl| ) and (22), the esti- 



mated error is smaller than 0{NApAt\\H\\/h)+0{At^\\H\\ max{\\H±\\} /h'^) 
+C'(iVApAt||max{||i/±||}/;i) +0{At^[[H+\\\\H_\\/h^). We also have ap- 
proximated the q-number denominator of the Gaussian operation elements 
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by a c- number, cp. (3C), which leads to an error of order C'(ApAt7). We 
further expanded the operation in powers of At up to the first order, cf. 
( p2[ ) and (Hil), which results in errors of order 0{At^\\H + HAvf/fi^) and 



0{At'^^\\H + HavW/^)- the continuum limit the errors all vanish but 
they can play an important role for discrete sequences of measurements if 
Ap ^ 1 is not fulfilled. Appendix C contains a calculation of the state 
change up to higher orders. This more accurate calculation confirms the 
order of the errors estimated here. In case of the neglected commutators in 
(^T|) and (|23|) the error turns out to be actually smaller (cp. (|68|)). 



5 Continuum limit 

In Section 4 we have worked out the operation given by the discrete state 
transformation between an initial state p{t) and the final state p{t + At) 
after a N-series of instantaneous interactions of a qubit with an environ- 
ment, which are of the type of an unsharp measurement. It can be applied 
to a truly sequential measurement by dividing the sequence of elementary 
measurements into a succession of N-series. Since the r.h.s. of the equation 
( ^SD is proportional to At, the given approximation is not sensitive to the 
division as long as N is large. The rate of the change of p is invariant under 
this the division. This discrete-time analysis is the most natural approach 
to sequential measurements. For example see [^. Eqn. (^) above reveals 
the underlying physics as represented by the decoherence rate 7 and the 
measurement induced unitary development given by H^v- 

There are elaborated schemes for the treatment of permanently open 
quantum systems by continuous-time descriptions. Master equations are 
an example. One may profit from these schemes as approximations in the 
sequentially open case too, if a physically reasonable continuum limit r ^ 
is carried out. The corresponding demand for such a limit is that the physical 
characteristics of the sequential situation have to be taken over. We proceed 
as follows: 

The quantity po is the mean probability to obtain the measurement result 
-|-. We leave the value of po unchanged in the continuum limit. In order 
to not change the decoherence behavior in the continuum limit we secondly 
demand for the decoherence rate: 

(Apf 

lim - — ; — = 7 = const . (49) 

T^o ipo{l - pq)t 

The smaller Ap the weaker the single measurement. With r ^ and the 
strength Ap of the single measurement unchanged, a Zeno effect would be 
obtained. This is prevented by appropriately diminishing the strength Ap 
of the measurement according to (p9|). This demand can also be found in 
the literature |T^ . 
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If in a given sequential physical situation the H^v is non- vanishing, then 
the total Hamiltonian dynamics is according to (^) governed by the Hamil- 
tonian H + Hj^y. We want to keep this dynamics in the continuum limit 
on physical grounds and demand therefore that Hav remains unchanged. 
Performing the limit r ^ as specified above results in the master equation 

P = -^[H + Hav,p]-IW,,Wz,p]], (50) 

which describes approximately the discontinuous situation in the noisy chan- 
nel characterized above. 

We note that the master equation ( |50|) could have been obtained in the 
limit ( ^9|) directly from the elementary measurements. The intervening for- 
mulation in terms of N-series will be exploited by the forthcoming equations 
of selective evolution. 



6 Selective evolution 

In Sec. 3 we calculated the Gaussian form ( p^ of effective operation ( |35| ) 
valid for a N-series. In Sec. 5 we derived the master equation (|50|) valid 
exactly in the continuous limit (p9D. As a matter of fact, the master equation 
describes the non-selective evolution. Selective evolution is, on the contrary, 
conditioned on the random measurement results (readout) and described by 
stochastic equations. In our case, the readout is s. It is the continuously 
measured unsharp value of the observable az obtained in the N-series in the 
limit (H). 

The theory of the selective evolution has been available since long ago 
15|. From the Gaussian operations (p4[) in the limit At — > 0, it has been 



proved that the selective evolution of the quantum state, conditioned on the 
measurement result s, satisfies the conditional master equation: 

P = -^[H + Hav,p]- lWz,Wz,p]]+w^{az - {az),p} (51) 

The function 'w{t) is the standard white-noise and the equation should be 
understood in the Ito-stochastic sense. The state evolution couples to the 
readout s by: 

s = (a,) + -^w (52) 

Obviously the stochastic mean of the conditional master equation ( [5l| ) re- 
duces to the unconditional master equation ([50|) as it should. Of course 
eqn. (51) applies to pure initial states as well. Then the pure state property 



p is preserved. The derivation may completely be identical to that in 



Ref. [15|. In the continuum limit (^) the value of Ap must vanish and the 



feed-back Ug is thus deterministic, given by Hav alone. 
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The above equations of selective evolution are exact in the following 
sense. Elementary operations are being applied with frequency growing to 
infinity and strength decreasing to zero as given by (^), i.e. at fixed 7. We 
read out the rate N^/N averaged over time At which should go to zero in 
such a way that N = At/r still goes to infinity. The elementary time r 
goes "faster" to zero than the time At to calculate the rate N^/N. The 
calculated current rate N^/N is related to s by (^): 



The continuous limit ( |49D of s exists. As follows from (52), it is centered 
around a state dependent part (cr^) and superposed by the white-noise of 
constant intensity I/7. 

Appendix A 

In the appendices we sketch the calculation of the change of state in the 
non-selective regime including all terms up to order 0{AtAp'^) and O(At^), 
where At occurs in products with either H/fi, H±/fi or 7. 

We start with the exact operation element for a N-series with unitary 
development U between consecutive measurements. 

■.= Um^ ...U^,\M^^ \ ...\M^,\U^ + Ri + R2 , (54) 

where Ri and R2 are the terms which arise from commuting out the evolu- 
tion operators U and the feedback operators U± respectively, cp. equa- 



tions (|2l| ) and (23). The relation of Ri and R2 to Ci and C2 and an 
estimation of the order of magnitude of Ci and C2 is described below. 
Since the commutators Km^ := [U , M^J occurring in Ri are of order 
0{ApT\\H\\/ti) + 0{t^\\H\\ max{||i^±||}/?i^), we can neglect terms contain- 
ing products of two such commutators. 

+ ... +(iV-l)M„^...K„,C/^-i (55) 



N-l 

n 



N-l 



ApT[H,a,]Um^...Um, \MmJ:.\MmAU 

^2 N-l 
n=l 



N-l 



where a+ := — 4y/po and a_ := 4-^1 — pQ. Please note, that in ( |55D behind 
the sum sign the products of |M|'s and C/'s are meant to not contain Mmpf_„ 
and Umj^_„, except if explicitly mentioned. In the last two lines we have 
commuted out the feedback U±. The resulting error is of higher order and 
can be neglected. 
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Also in R2 we only take into account the terms containing one commu- 
tator Kk^i := [Mm,. , Umi ], which is of order 0{IS.pT ma,-x{\\H±\\} /%). 



R2 = Umj,.:Um,K2,l\MmJ...\Mm,\\MrnAU^ + - (56) 



N 



X \M.^J ...|M^,^J|M^,_J ...|M^jC/^-i 

X ;7™,_^ ...Um,\M.raJ ...\Mm,\U^-'^ ■ 

Let us shortly motivate the estimation of the order of magnitude of Ci and 



C2 given in (|2|) and (||). First we observe that Ri = Mm„ ...Mm^U^Ci 
and R2 = uf+u!^''^+\M+\^+\M_\^-^+U^C2. A moment's thought shows 
that the order of magnitude of the summands contained in Q is equal to the 
order of the commutators and K^^i in Ri. Since there are approximately 
iV^ such summands in Ci, the norm of Ci can be estimated to be less or equal 
to A^^ times the order of the commutators in Ri which leads to the claims 
(I2I) and (|24 



Appendix B 

The state change due to a N-series in the non-selective regime reads 

Pit + At) = J2^(^n,)P^U) 

J2 Um^ ...U.raAMmJ-\MmAU(At) X 
mi ... mjv 

xpU+{At)\M,mA■■■\Mrr^^\U+^...U+^ 

+Ri + R2 + 0{Rl) + 0{Rl) + 0{RiR2) (57) 

with 

Ri:= {RipU+{At)\MmA...\MmAU^i-U^M+h-c] i = l,2. 

mi ... mjv 

(58) 

In order to carry out the summation in (^7|) its terms can be expressed 
by means of binomial distributions. The latter can be approximated by 
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integrals over Gaussians. For the first term in (|5^) this recipe has already 
been demonstrated in section 3 and 4. Ri can be written in binomial form 
by observing that: 

mi ...ruM n=l 

r+ 



Xp \MmA-\MmJUZ,-U, 



rriN 



X 



\ m / N+=0 

x\MXMM^\'^-''+-\\M+f+\M.f-^+-\ut)^+{U+)''-''+-\ (59) 

where we have again used the notation that operators with upper script L 
and R act from the left and from the right respectively. A similar formula 
is obtained when instead of |Mmjv_„| in the first line UmM~n+i missing. 
Then only has to be replaced by |Mm|^. 
i?2 can be simplified employing 

N k-1 

X/ X] ^rnu,miUmN ■■■ Umk+iUnik-i ■■■ Umi X (60) 
mi ...m]v k=2 1=1 

x\Mmr,\ - \Mmi+A\Mm,_A- l^mi | P ... |M^^|C/+^ ... [7+^ 

\rn,rh=+,~ 

N+=0 \ + / 

X p |M+ 1 ^+ 1 M„ I ^-^+-2 ([/+ )^+ ([/+)^-^+-2 , 

Formulae ( ^9|) and (^) neglect commutators between the operators they 
contain. In our case corrections containing these commutators would be of 
higher order and therefore too small. 

Applying formulae (|5^), (|60|) to Ri, R2 respectively and expressing 
A'^^ /N in terms of variable s according to ( |3ll ) we obtain: 



Ri = -^^^^\[H,a^]Y.^D{p)\M^\U;^ + h.c.\ (61) 



At2 



J2{[H, Hrn]\Mra\D{p)\Mm\U+ + h.C.} 



where Dip) = J^^MspM^ds with Ms as given by (|3^) with At replaced 
by (N — l)r. R2 now reeds 

i?2 = - '^2h^^ ^ | ^'^'^^'^V ^|M^|£)(p)|Mrf,|[/t|M^|^+ + /i.c.| 
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E {[Hm, H^^A\Mm\\M^\D{p)\M^\UX\Mm\U+ + h.c.} (62) 

In i?2, when inserting Mg from (|35| ) in D{p), At has to be replaced by 
(A^ — 2)r. The second sum in ( |6^ ) vanishes since the summand with m = 
+ , m = — and the summand with m = — , rh = + add to zero. Inserting 
the lowest order of D{p) namely D{p) p, it is easy to see that Ri and 
i?2 contribute the terms in the sixth and the seventh line of equation ( |6^ ) 
respectively to the change of state. 



Appendix C 

Having calculated i?i and R2, we want to sketch how to process the main 
contribution to the state change, which is represented by the first term in 
equation (|57[). As mentioned above this part of the operation can be written 
by means of a binomial distribution and then be expressed with operation 
elements whose modulus |M(A^+,A^)| are the square root of Gaussians, cf. 
equation (^). In contrast to section 3 and 4, we now take into account the 
full q-number denominators of the modulus |M(A^+, A^)| in (29). Expressing 



the operation elements in terms of variable s (cp. (pl|)) we obtain for their 
modulus 

|MJ = ^^^=exp|-7^^^-^Atl, (63) 
' ^27r/(7At) ^1^4 j' ^ ^ 

with 

(64) 

Expanding the unitary part of the operation up to order At^ leads to the 
state change (without the contribution from R\ and R'l) 

iAt At^ 
p{t + At) = D{p)- — [H + HAv,D{p)]-^{{H + HAvf,p} 

At^ 

+-^{H + Hav)p{H + Hav) (65) 
n 

with 

D(p) = y_%iip|^(A/;L-AHR)sApAt|exp|-^(aL-,)'At| X 

We note that in eqn. (35) H is meant to act in operator products directly 
on rho. This is due to the order of operators in the operation elements (cp. 
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(|36D). The integral D{p) has a closed form solution which can be expanded 



in powers of At and Ap. 7 without hat is given by (33). 



1 - At 



2(.L_,R)2^1_lfMP0-l/2) 



(67) 



2 V PoPo 
X (1 - Vp^iPoPo + 3 - 2-1/4))) 

+^ (aH^ - AH^Y + ^^{AH^ - AH^Xa}' + af 

+At'^^ {a} - aff + 0{ApAt'^) + ©(Ap^At^) + ©(At^) I p . 
32 I 

Collecting all terms up to order C'(At^) and O(AtAp^) we obtain the 
following difference equation: 



Ap = At\^--[H + HAV, p] 

_7 ^Ap(po-l/2) 



2 V PoPo 



(68) 



(1 - Vp^{P0P~0 + 3 - 2-1/4)) j [^^ ^ [^^ 



^[AH, [AH,p]] -!^[AF,K,p}] 



87?!^ 



4/i 



+At^\^-^{[a,,[a,, [H,p]]] + [Hav, , , p]]]) 

-^{{H + Hav?,p} +^{H + Hav)p{H + Hav) + ^Wz,Wz,p]] 



iAt^^ SiAtAp 



+ 



(^l^[H, a,] pa, + h.c.^ -^[[H, Hav_ 



\ 2h shpopQ 

+^^[(^z , PoH+ + pqH^] paz +h.c.|j 

0{NApAt'^) + 0{ApAt^) + 0{Ap'^At^) + 0{At^) 

with pq := 1 — pq. In the order terms At occurs in products with one of the 
three: H/h, H±/?i or 7. 
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